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ABSTRACT 

A new classification scheme of the IRAS LRS carbon stars is presented. It comprises the separation of 718 
probable carbon stars into 12 distinct self-similar spectral groupings. Continuum temperatures are assigned 
and range from 470 to 5000 K. Three distinct dust species are identifiable: SiC, a:C— -H, and MgS. In addi- 
tion to the narrow 1 1 + pm emission feature that is commonly attributed to SiC, a broad 1 1 + pm emission 
feature, that is correlated with the 8.5 and 7.7 pm features, is found and attributed to a:C— H. SiC and 
a:C— H band strengths are found to correlate with the temperature progression among the Classes. We find a 
spectral sequence of Classes that reflects the carbon star evolutionary sequence of spectral types, or alterna- 
tively developmental sequences of grain condensation in carbon-rich circumstellar shells. If decreasing tem- 
perature corresponds to increasing evolution, then decreasing temperature corresponds to increasing C/O 
resulting in increasing amounts of carbon rich dust, namely a:C““H. If decreasing the temperature corre- 
sponds to a grain condensation sequence, then heterogeneous, or induced nucleation scenarios are supported. 
SiC grains precede a:C— H and form the nuclei for the condensation of the latter material. At still lower 
temperatures, MgS appears to be quite prevalent. No 11.3 p m PAH features are identified in any of the 718 
carbon stars. However, one of the coldest objects, IRAS 15048 — 5702, and a few others, displays an 11.9 pm 
emission feature characteristic of laboratory samples of coronene. That feature corresponds to the C— H out 
of plane deformation mode of aromatic hydrocarbon. This band indicates the presence of unsaturated, sp 3 , 
hydrocarbon bonds that may subsequently evolve into saturated bonds, sp 2 , if, and when, the star enters the 
planetary nebulae phase of stellar evolution. The effusion of hydrogen from the hydrocarbon grain results in 
the evolution in wavelength of this 11.9 pm emission feature to the 11.3 pm feature. 

Subject headings: circumstellar matter — infrared: stars — molecular processes — stars: carbon — 
stars: fundamental parameters 

1. INTRODUCTION 

Historically, silicon carbide (SiC) was the first dust species 
identified from its infrared spectroscopic features in carbon 
stars (Treffers & Cohen 1974). The 1 1 pm band profile has been 
observed from ground based observatories in a large number 
of carbon stars by Forrest, Gillett, & Stein (1975) and by 
Cohen (1979). Whereas Forrest et al. (1975) found a fairly 
uniform band profile among the visually bright carbon stars, 

Cohen (1980) observed subtle differences in cold dust 
enshrouded carbon stars, those he termed extreme carbon 
stars. In the bright visible carbon star Y CVn, Goebel et al. 

(1980) observed a somewhat different SiC band profile than 
normally found in the other visible carbon stars that was com- 
parable to that expected from the crystalline forms of SiC pro- 
duced in the laboratory by Dorschner, Friedemann, & Gurtler 
(1977a, b). Until the IRAS Atlas became available, there were 
available for study relatively few spectra of the 1 1 pm SiC 
band. With the advent of the LRS, there are hundreds of 
spectra with good signal-to-noise ratio and spectral coverage 
in the important 15 pm region that is blocked from ground and 
airbomed observation. The characteristics of the band profiles 
have been the subject for discussion by Papoular (1988), Little- 
Marenin et al. (1987), and Skinner & Whitmore (1988a, b). 
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SiC has been incorporated into models and is currently 
studied in the laboratory as a carbon star circumstellar 
material. SiC particles have been studied extensively in the 
laboratory by Dorschner et al. (1977a, b), Friedemann et al. 
(1981), Borghesi et al. (1983), Borghesi et al. (1985), Pegourie 
(1988), and Lewis, Amari, & Anders (1990) with the astro- 
physical application in mind. Radiative transfer models for 
stellar carbon rich circumstellar shells concentrated mainly on 
the carbon dust component because of the significant absorp- 
tion coefficient it displays at nearly all infrared wavelengths, 
and treated SiC as an impurity by mass and a minor contrib- 
utor of opacity, and therefore, to the balance of energy in the 
radiation transfer equation (Martin & Rogers 1987). In most 
modeling of carbon stars, the SiC feature and component by 
mass is relegated to insignificance (e.g., Rowan-Robinson & 
Harris 1983; Le Bertre 1988). To the contrary, the models of 
Chan & Kwok (1990) presume only an empirical SiC opacity 
without carbon dust that can model the radiative transfer suc- 
cessfully in optically thin and thick shells. The peak in the 
spectral energy distribution of some planetary nebulae occurs 
in the 8 to 14 pm region, so that SiC emission feature can 
represent a significant fraction of the total radiation in the 
spectrum. See Hoare (1990) for a recent application of dust 
models that examine the dust-to-gas ratio in carbon-rich plan- 
etary nebulae. 

The observational studies of the SiC component in meteor- 
itic specimens by Tang & Anders (1988), Zinner, Tang, & 
Anders (1989), and Amari et al. (1992) indicate the presence of 
SiC outside the confines of a carbon-rich circumstellar 
environment and furthermore are of importance to our under- 
standing of the early solar nebula. The larger meteoritic SiC 
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grains contain rare gas isotropes consistent with those that are 
expected to be present in the asymtotic giant branch (AGB) 
He-burning shell and the envelope (Lewis et al. 1990; Gallino 
et al. 1994; Hoppe et al. 1994). The smaller meteoritic SiC 
particles also appear to be possible carbon star ejecta with 
s-process enhancements of isotopes (Brown & Clayton 1992). 
Another SiC component with unusual isotopic composition, 
grains “X,” may be formed in the explosive ejecta of super- 
novae (Amari et al. 1992; Hoppe et al. 1994). All the SiC grains 
examined by Hoppe et al. (1994) are cubic (0-SiC) based on 
X-ray diffraction, yet amazingly many different grain forms are 
displayed. 

Historically, the question of the presence of SiC in the inter- 
stellar medium has been of considerable interest, and its puz- 
zling absence in the interstellar extinction curve is the subject 
of a recent discussion by Whittet, Duley, & Martin (1990). 
They set an upper limit on the mass fraction of SiC in inter- 
stellar dust of at most a few percent, arguing that it is oxidized 
in the interstellar diffuse cloud environment on a timescale of 
~5 x 10 7 years. 

It is the purpose of this work to demonstrate that there is a 
progression of 11 /mi emission features in carbon rich stars 
that indicate the presence of material components other than 
SiC. Spectral classes are defined by the use of automatic classi- 
fication procedures and can be understood in terms of multiple 
grain types. It is possible to find reasonable spectral identifica- 
tions for the emission features that are consistent with physical 
conditions in and evolutionary history of carbon stars. 

2. A NEW CLASSIFICATION OF LRS CARBON STARS 

Previously, we used the artificial intelligence classification 
technique AutoClass II developed by Cheeseman et al. (1988) 
to automatically classify all the spectra in the Atlas of Low- 
Resolution IRAS Spectra (1986, hereafter the IRAS LRS 
Atlas). This is in contrast to the human expert technique used 
by the authors of the IRAS Atlas that was based on the 
obvious concept of spectral similarity. The study upon which 
we report here is an extension of the previous Cheeseman et al. 
(1988) work and focuses only on the groups identified as 
carbon star Classes. It includes 718 objects which we will call 
probable carbon stars. The result of this effort is summarized in 
Table 1 ; we find a total of 12 separate classes of carbon stars in 


the LRS. AutoClass II has produced 12 carbon star Classes 
from a sample of 718. A total of 10 or 1 1 classes are expected by 
application of Sturges’s rule for classifying a series of N items 
(Newman, Haynes, & Terzian 1994). AutoClass II appears to 
be finding statistically significant differences in the population, 
not arbitrary groupings. 

The objects are drawn from the MetaClasses al-a6, SI, 220, 
and 230:1 and 230:5 of Cheeseman et al. (1988). Also included 
were all those carbon stars identified in the catalog association 
field (catalog 17; Stephenson 1973) of the IRAS LRS Atlas 
distribution tape even if not included in the carbon star Meta- 
Classes, e.g., BM Gem, which is in a silicate emission Class (11). 
It would be silly to intentionally exclude visibly identified 
carbon stars, for that is totally inconsistent with the Bayesian 
principles of incorporating all prior probabilities. In retrospect, 
we have missed a small number of stars which did not satisfy 
these criteria like 08445 — 3932 = HD 75021, since we assumed 
all visibly bright identified carbon stars were included in the 
CCCS = Catalog 17. Hence, we are not claiming absolute 
completeness for the list of carbon stars presented here. Not 
included were MetaClasses y and (, among which are found 
some protoplanetaries and many planetaries, since we did not 
wish to expand the working definition of carbon stars to 
include both those groupings. Moreover, their spectra are suffi- 
ciently similar to those of silicate absorption shells so that the 
LRS instruments give little distinction between the two groups. 
The reference work by Cheeseman et al. (1989) provides a 
separate listing of planetary nebulae and protoplanetary 
objects from the LRS that are primarily resident in the y 
MetaClass. 

Considering only the MetaClass members, we find 161 out 
of 176 stars found in the list of Little-Marenin et al. (1987) of 
new carbon stars and 57 out of 65 of the carbon stars listed by 
her as well known. When catalog 17 association is included, all 
their previously determined carbon stars in the LRS Atlas are 
recovered. In addition, there are objects in our list that are new 
associations with carbon stars, determined solely upon the 
basis of the AutoClass search reported here. Most of these 225 
probable carbon stars are of low signal-to-noise ratio. Their 
validity needs to be checked at near infrared wavelengths by 
spectroscopic means. We believe that the vast majority of 
carbon stars have been included in this study, and that the few 


TABLE 1 


Summary of AutoClass II Carbon Star Characteristics 


Class 

Population 

% 

(K) 

Dust/Gas Specie s 

LRS Classes 
(Dominant) 

Galactic 

Latitude 

Dispersion 

Mag 

[12] 


FIR Colors 


[12]— [25] 

[25]— [60] 

[60]— [ 100] 

0 

109 

800 

SiC + oc:C — H 

4n 

15° 234 

-1.053 

0.462 

0.229 

|Kta 

1 

103 

600 

a:C-H + SiC + MgS 

4n 

6.344 

-1.119 

0.867 

0.256 


2 

74 

5000 

Unknown 

On, In 

14.288 

1.070 

0.528 

0.328 

2.511' 

3 

74 

2000 

Gas + dust 

4n, In 

20.679 

-0.428 

0.262 

0.531 


4 

73 

2000 

SiC 

4n, In 

9.365 

0.324 

0.492 

0.342 

1.113 

5 

60 

470 

a:C — H + MgS 

4n, In 

7.739 

-1.016 

1.222 

0.395 


6 

58 

5000 

Gas + dust 

In, On, 4n 

11.638 

0.764 

0.344 

0.456 

1.735 

7 

46 

5000 

Gas + dust 

On, In, 4n 

7.929 

1.175 

0.419 

0.404 

2.140 b 

8 

39 

900 

SiC + a:C— H 

In, 4n 

4.256 

0.790 

0.653 

0.374 


9 

39 

600 

ot:C — H 

In, 4n 

8.547 

0.554 

0.570 

0.533 

0.9 60 c 

10 

29 


Unknown 

On, 2n 

12.071 

1.373 

0.731 

0.748 

1.861** 

11 

6 

600 

Gas/silicates 

2n, In, 4n 

8.119 

0.318 

0.839 

0.022 

0.181 


1 Only two stars with reliable [60]— [100]. Must be considered an upper limit. 
b Only three stars with reliable [60]— [100]. Must be considered an upper limit. 
c Only four stars with reliable [60]-[100]. Must be considered an upper limit. 
d Only one star with reliable [60] [100]. Must be considered an upper limit. 
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that have been overlooked will not greatly affect the conclu- 
sions based upon looking at the groups, rather than individ- 
uals. 

The parameters of the 12 carbon star Classes are sum- 
marized in Table 1. “Class” is the AutoClass Class number; 
“Population” is the Class membership count; T c is the esti- 
mated Class continuum temperature; “ Dust/Gas Species” lists 
the dust species assigned here; “LRS Classes” lists the domi- 
nant LRS Classes present in an AutoClass Class; “Galactic 
Latitude Dispersion ” is the standard deviation in the distribu- 
tion of galactic latitude of members of the Class; “ mag[12] ” is 
the IRAS band 1 flux magnitude average of the Class as 
defined in Walker & Cohen (1988) and Cheeseman et al. 
(1989); and “ FIR Colors ” are the IRAS band 2, 3, and 4 color 
differences of the Class average in magnitudes as defined in 
Cheeseman et al. (1989). 

The temperatures quoted in Table 1 are derived by compari- 
son of the flux curves with single temperature blackbody func- 
tions, as distinguished from a more complex radiative transfer 
calculation involving dust opacities. The fitting criteria are 
unique to each Class taking into account the different charac- 
teristic spectral features of each Class. Temperature assignment 
represents post-AutoClass analysis. Of course noisy Classes 
have less reliable temperature estimates of the continuum, T ci 
than quiet Classes. The absolute temperature scale is some- 
what suspect because of the photometric calibration proce- 
dures in the original IRAS LRS Atlas which used a Tau as the 
primary standard and assigned to it a temperature of 10,000 K, 
as the accepted temperature of a Tau is more like 4500 K. We 
have not attempted to correct for the temperature calibration 
error present in the data from the Atlas in order to maintain 
consistency with the data set presented there. Due to the pres- 
ence of SiO absorption in the 8 ji m region, the red end of each 
and every spectrum is about 10% higher in flux than it should 
be. Approximate adjustment can be made by using the ratio of 
10,000 and 4500 K blackbodies in the appropriate fashion. 
Recalibrations are discussed by Wooden et al. (1993). The 
more exact correction can be made using the recalibration 
table of Cohen et al. (1992). Our choice to maintain consistency 
with the LRS atlas is not without penalty. We may have missed 
identifying those stars that have intrinsic SiO absorption at 8 
pan. SiO is predicted to be abundant in certain elemental abun- 
dance circumstances in carbon stars (Tsuji 1964, 1973), and has 
been identified at 4 /mi in the spectrum of TX Psc by Bregman, 
Goebel, & Strecker (1978). 

In Figure 1, we compare the spectra of Classes 0, 1,3, 4, and 
5 with the laboratory spectra of SiC and a:C— H. The Classes 
have been arranged from top to bottom in order of increasing 
8.5 /im emission strength and decreasing 11-}- /mi emission 
strength. This sequence corresponds to a continuum tem- 
perature, T c , sequence decreasing with increasing 8.5 /an emis- 
sion excess, with the exception of Class 3. Analysis of Class 3 
shows that it is complicated by the presence of molecular 
absorption bands (Goebel et al. 1988), so we can understand its 
deviation from the T c sequence on this basis. We consider there 
to be a good T c sequence for the Classes which are apparently 
dominated by dust, i.e., 0, 1, 4, and 5. 

From the comparison of Class spectra with material emit- 
tances that is displayed in Figures 2 and 3, we can conclude 
that SiC predominates in the higher temperature spectra, while 
cr.C — H dominates at the lower temperatures. At intermediate 
temperatures, 500 K < T c < 2000 K, and excluding Class 3 for 
the stated reasons, there appears to be a mixture of the fea- 



Fig. I. — The average spectra of Classes 0, 1,2, 3, 4, and 5, the most popu- 
lous carbon star Classes, normalized to their respective continuum tem- 
peratures. Also displayed are the arbitrarily normalized absorption coefficients 
of /?-SiC and a:C— H for comparison purposes. 


tures. The mixture causes the peak wavelength and width of 
the 1 1 -b /mi feature to shift from Class to Class. We conclude 
that the shift cannot be due to temperature effects on the SiC 
emission band, because the predicted broadening of the band 
profile increases with increasing temperature. This is in oppo- 
sition to the observed trend. In Table 1, we have listed the dust 
components that appear to be correlated with each of the 12 
Classes. Some of the Classes, like 8 and 9, are low S/N versions 
of Classes 0, 1,4, and 5. 



Fig. 2. — The average normalized emission excess spectrum of carbon star 
Class 4 compared directly with the absorption coefficient of fi-SiC normalized 
to the same 1 1 4- /im peak height as Class 4. SW and LW refer to the short and 
long wavelength data sets of the IRAS LRS spectrometer. 
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Fig. 3. — The average normalized emission excess spectrum of carbon star 
Class 5 {solid and dashed curves ) directly compared with the absorption coeffi- 
cient of a:C— H ( dotted curve ) normalized to the same 1 1 + jim peak height as 
Class 5. 

2.1. Other LRS Classifications of Carbon Stars 

Little-Marenin et al. (1987) have classified 176 LRS spectra 
as carbon stars using the 1 1.2 /mi SiC emission as the selection 
feature. They compared an appropriate blackbody function to 
each spectrum and examined the remaining excesses for the 
characteristic 11.2 /im emission profile. Our examination has 
recovered all of the bona fide entries in their Tables 1 and 2 of 
known and previously unknown carbon stars. Outside the 
scope of the present discussion are the puzzling few carbon 
stars with silicate-like emission features. Both Little-Marenin 
(1986) and Willems & de Jong (1986a, b) have discussed the 
discovery and implications of these spectra. AutoClass II also 
recognizes these spectra and places them in a separate Class 1 1 
of Table 1. 

Willems (1987) considered the 304 carbon star spectra listed 
in the LRS Atlas. He provided an SiC index for those with an 
emission feature peaking between 11.2 and 1 1.6 /mi. He recog- 
nized spectra with featureless dust emission outside the SiC 
feature attributed to amorphous carbon dust. Other spectra 
exhibited absorption bands due to molecular species, primarily 
HCN and C 2 H 2 . Those with near-infrared color temperatures, 
7 J r > 2000 K, have optically thin dust emission shells in the 
continuum and frequently display the SiC feature as well as 
molecular absorptions; while those with 7J R < 2000 K have 
optically thick dust emission shells in the continuum and 
exhibit emission from dust, either from SiC or from amorphous 
carbon, or both, along with molecular features. Since the 
maximum SiC index was 0.24, the interpretation was given that 
only a limited fraction of the carbon was condensed in the form 
of SiC, of order l%-2% of the dust particles. The variation in 
shape of the SiC feature was considered to be consistent with 
the well-known sensitivity of the feature to particle shape, as 
was discussed by Treffers & Cohen (1974). The limit of the 
classification effort to optically thin and thick continuum 
shells, largely left unaddressed the issues to be discussed here, 
namely the progression of emission features within the LRS 
sample. 


Popular (1988) has classified the infrared spectra of carbon 
stars in the LRS following a different procedure than ours. He 
observed three different types of 11.3 /im features differing in 
width and by the occurrence of secondary features at ~8.6, 
~11.7, and 12.8 /im. The different SiC features could be 
matched by pure or polluted a-SiC from the laboratory studies 
of Bohrgesi et al. (1985). The classification procedure was that 
of Baron et al. (1987). First an underlying continuum was esti- 
mated based on a power law connecting two spectral points in 
the short-wavelength portion of the spectrum as chosen by the 
observer. The resulting spectral excess and feature shape were 
defined mathematically. Then the individual spectra were 
ranked according to excess and shape of the SiC feature, 
nominally at 11.3 /im. The resultant spectra displayed a wide 
range of behaviors in the long wavelength region that were 
primarily due to the estimation of the continuum by a power 
law rather than a blackbody. However, in the short- 
wavelength region, the feature shapes have much in common 
with the feature shapes in the AutoClass II Classes of the 
present work. For example, his subset SiC(a) is similar to our 
Class 0 and subset SiC(b) is similar to Class 4. There exist 
significant differences as well. Subset SiC(c) has no counterpart 
in our study. The closest we come to such a spectral shape is 
Class 1 1, the Class which contains the few stars with silicate 
type emission. While we also use an estimate of the continuum 
and remove it to study the spectral features, we chose black- 
body functions because of the significant spectral curvature in 
the 8 to 24 /xm wavelength range of the functions in the lowest 
range of temperatures appropriate to these stars. 

The classification by AutoClass II was done in flux versus 
wavelength space, not blackbody normalized excess versus 
wavelength space, so that the temperature of the objects would 
be a consideration in the procedure. Papoular’s procedure 
removes the influence of temperature prior to classification. 
Therefore any influence of temperature on the spectral features 
is lost. Additionally, the starting set of objects in the two 
studies are significantly different. Our set was selected pri- 
marily based upon the previous AutoClass II examination of 
the entire LRS Atlas, whereas Papolar’s set was selected by the 
eye of an expert individual, rather than a mathematically rigor- 
ous algorithm. For these reasons, our results are significantly 
different. 

Chan & Kwok (1988, 1990) studied 356 carbon stars from 
the General Catalog of Cool Carbon Stars (Stephenson 1973) 
by examining the IRAS far-infrared photometry and the LRS 
spectrophotometry of the 538 objects listed in LRS Class 4n. 
The sample was paired down to 356 for various reasons. The 
extent of their classification study was to redefine the spectral 
index with the result that none of the SiC emitting spectra had 
indicies exceeded their index 45. They studied the differences 
displayed in the characteristics of the SiC emission feature in 
the visually bright and infrared bright carbon stars. The differ- 
ences in the SiC spectra of the two groups were interpreted to 
express an evolution from a-SiC to /?-SiC parallel to the evolu- 
tion from visually bright to infrared bright carbon star. Our 
approach to classification and spectral analysis differ signifi- 
cantly from their work, and we come to distinctly different 
conclusions regarding the nature of the dust emitting 
materials. 

Volk & Cohen (1989) examined the LRS database from 
which the LRS Atlas was extracted, and found 365 high- 
quality spectra that for one reason or another were not 
included in the IRAS LRS Atlas. Of the reported sources, 28 
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were determined to be SiC emitters and some dozen others 
possibly carbon stars. The large number of potentially carbon- 
rich planetary nebulae that they found are not considered here. 
Of the nearly 40 carbon stars outside the LRS Atlas, some are 
well-known visibly bright carbon stars, e.g., TX Psc, V CrB, 
S Cep, and V Hya. They have not been included in our analysis 
of the Atlas, because we examined the IRAS distribution 
product, or tapes, not the hours-confirmed-sources, or raw 
data distribution on CD-ROMs. Several of these extra sources 
are included in the study by Skinner & Whitmore (1988a, b) of 
the SiC feature and its relation to mass-loss rates. 

3. DISCUSSION 

We divide the discussion into three sections. The first section 
discusses the similarities shown by the LRS spectra with labor- 
atory spectra of relevant materials and the dissimilarities of 
other LRS emission features. The second section discusses the 
evolutionary theories proposed here and elsewhere. The third 
section deals with the implications of the present interpretation 
for dust condensation theories. 

3.1. Laboratory and Astronomical Spectra 

In Figure 2, we present a direct comparison of the SiC emis- 
sion with that of Class 4. The a-SiC emission is from the work 
by Pegourie (1988). The 300 K temperature of the SiC band 
produces an understandably narrower emission profile than 
the 2000 K continuum of Class 4 might imply. While it is not 
obvious what temperature to assign to the emission of Class 4’s 
SiC profile, it appears to be warmer than 300 K, and, therefore, 
broader than the laboratory sample, specifically on the long 
wavelength side, rather than the short wavelength side of the 
feature. The sharp cutoff on the short wavelength side of the 
feature disallows a strong temperature influence on that part of 
the feature and accounts for the similarity of that profile. The 
slower cutoff on the long wavelength side allows the tem- 
perature to influence the feature width. It would be interesting 
to model the shape and temperature of the SiC grains to deter- 
mine the effect of these parameters on the profile of the band. 

The a:C—H absorption spectrum in Figures 1 and 3 is from 
Dischler et al. (1983a, b). It compares most favorably with the 
spectrum of Class 5, the Class with the lowest continuum tem- 
perature assignment, 470 K. The direct comparison of Class 5 
with a:C — H is presented in Figure 3. The reported laboratory 
absorption shown in Figure 1 displays stronger absorption at 
7.8 jun than at 1 1 4- pm. Normalizing the laboratory band- 
strengths to a 300 K excitation temperature gives nearly equal 
bandstrengths. In Figure 3, the temperature normalized bands 
of the laboratory material and the temperature normalized 
bands of the stellar spectrum compare very well in central 
wavelength and symmetry. The central wavelength of a:C—H 
is 11.55 pm, and that of the feature in Class 5 is 11.67 pm. 
Minor differences are apparent in the width of the feature at 
11 -I- pm and the shape of the base of the feature. In Figure 1, 
the sharp turndown in LRS spectra short of 8 pm, that is an 
artifact of the LRS instrument, has been corrected using the 
recalibration of Cohen, Walker, & Witteborn (1992). Since the 
correction to spectra is greatest near 8 pm , and that region 
involves the 7.8 pm band of a:C— H, the corrected spectrum of 
Class 5 is used in Figure 3. 

As the Class temperature increases in Figure 1, so does the 
apparent strength of the 11+ pm emission relative to that at 8 
pm, with the exception of Class 3, about which more will be 
said later. The central wavelength of this feature also moves 


progressively to shorter wavelengths with increasing tem- 
perature. This can be interpreted in terms of a two-component 
model that mixes SiC and a:C— H, and favors SiC at higher 
temperatures and a:C~H at lower temperatures. The inter- 
pretation here is distinctly different from that arrived at by 
Papoular(1988). 

The apparent similarity of the two 1 1 + pm peaks for the 
two materials SiC and a:C— H leads to a natural explanation 
for the narrow and wide 1 1 + pm features observed by Cohen 
(1979) in carbon-rich dust shells. The narrow feature is 
attributable to SiC and the broad feature corresponds to 
a:C— H. SiC’s emission feature cuts off at a shorter wavelength 
on the long wavelength side of its emission feature than does 
that of a:C— H, see Figures 2 and 3. The alternative explana- 
tions for differences in the 1 1 + pm feature involve shape and 
size differences in SiC grains that influence the radiation cross 
section, or radiation induced damage effects on the crystalline 
structure. Because there are so many possible configurations of 
grain shape, size, and damage parameters, these explanations 
should produce a continuum of classes, whereas, in our 
analysis we observe a few distinct classes. 

The four Classes 0, 1, 4, and 5 represent 49% of all the 
AutoClass carbon stars, so we have an explanation for the 
characteristic dust emission profiles of a major fraction of all 
carbon stars in terms of two species. Classes 8 and 9, with an 
additional 1 1 % of the population, are of similar characteristics 
but with lower signal to noise ratio. Hence, the majority of 
carbon stars, at least 60%, can be described by a duality of dust 
species, a-SiC and a : C— H. 

A third type of 1 1.3 pm feature occurs only in the planetary 
nebulae and H n regions where there is significant ultraviolet 
radiation present. It is one of the well known and long studied 
UIR emission features. The 1 1.3 pm band is narrower than the 
SiC or oc:C— H emission bands. It has been attributed to poly- 
cyclic aromatic hydrocarbons (PAHs) (Leger & Puget 1984; 
Allamandola, Tielens, & Barker 1985). The LRS Class of these 
objects is 8n. Only two 8n objects, both members of Class 10, 
are found in the AutoClass carbon stars. The spectra of these 
two objects resemble hot continua with cosmic ray spikes 
rather than true 11.3 pm features. They may not be bona fide 
carbon stars. The lack of 8n Class objects in our sample arises 
from the preselection filtering of AutoClass MetaClasses done 
in the work of Cheeseman et al. (1989) that forms the basis of 
the present AutoClass carbon star set. The a MetaClass is 
predominantly populated by known carbon stars, i.e., those 
with catalog 17 (the CCCS) associations. These objects tend to 
have warm temperatures greater than 400 K and not the lower 
temperatures associated with the infrared spectra of planet- 
aries. Neither selection of objects specifically precluded LRS 8n 
membership, so apparently, the narrow 11.3 ^m feature is not 
associated with the circumstellar shells of carbon stars warmer 
than 400 K. 

It is noteworthy that hydrogen deficient R CrB type stars, 
which have enough ultraviolet radiation to be observable by 
WE, and are present in the AutoClass carbon stars, do not 
display the narrow 11.3 pm or the 8.5 pm emissions that are 
indicative of bounded hydrogen (Forrest 1974; Treffers et al. 
1976; Goebel et al. 1981; IRAS LRS Atlas). Therefore, the 
R CrB form of carbon dust must be more graphitic (Hartmann 
& Apruzese 1976) and deficient in hydrogen — and therefore, 
featureless in the infrared — than the hydrogenated tetra- 
hedrally coordinated, or diamond like, a:C— H found in the 
majority of carbon stars that display the 8.5 and 11+ pm 
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emission features. The graphitic form may be crystalline, or 
amorphous, or glassy, and must certainly be hydrogen poor. 
Without attached hydrogen atoms, there is no possibility of an 
11.3 /im emission, even with ultraviolet excitation. A lack of 
hydrogen in R CrB and other related hydrogen deficient stars 
is well known from other studies (see Warner 1967; Fadeyev 
1983; Walker 1986; Buss et al. 1993). 

The narrow 11.3 pm feature is currently associated with 
PAHs (Leger & Puget 1984; Allamandola et al. 1985). Our 
association of the broad 11 + /im band with a:C— H does not 
preclude the possibility that PAHs coexist with the dust, nor 
that a spectrum of particle sizes ranging from molecules to dust 
grains exists in the circumstellar shell. In the work of Goebel 
(1987) a mixed population of PAHs, HACs, and graphitic 
carbon were proposed to explain the infrared spectrum of 
planetary nebula NGC 7027. Subsequent studies by Bregman 
et al. (1989), Witteborn et al. (1989), and Buss et al. (1993) have 
added further support to the idea of coexisting populations. In 
NGC 7027 the narrow features were attributed to PAHs and 
excited by ultraviolet to a high temperature like 1000 K, while 
the broad underlying features (sometimes referred to as 
pedestals) were attributed to a:C— H and emitted at a moder- 
ate temperature like 300 K. The graphite component 
accounted for the far infrared continuum, which is featureless 
and radiates at 90 K with A’ 2 emissivity (Forrest, McCarthy, 
& Houck 1981). The ultraviolet radiation necessary to excite 
the narrow features like that at 1 1.3 /im is not present in the red 
giant phase of carbon stars. Even when there is a hot binary 
star present, as is thought to be the case for UV Aur, no narrow 
features are observed. None of the binary carbon stars in the 
list of Richer (1970) are in LRS Class 8n. Apparently, binary 
carbon stars do not have hot enough, or intense enough ultra- 
violet radiation fields to excite the features. The classical 
carbon stars do not display the narrow emission features. If the 
PAHs are abundant in these objects, they have to emit 
copiously at 11.3 jim, just like any other species must emit that 
has permitted vibrational bands in the infrared. Since the 
PAHs are essentially molecules, we can count the number of 
photons emitted by them at 11.3 jxm, and measure their popu- 
lation. The result is essentially nil. This observation certainly 
limits the abundance of PAH’s in the atmospheres or circum- 
stellar shells of classical carbon stars. 

On the other hand, one of the coldest objects of the present 
study (see Table 2), 15084-5702, Class 5, which is unas- 
sociated with any major catalog, shows the 11.85 /xm narrow- 
band characteristic of some laboratory PAH samples such as 
coronene (see Fig. 4). The 1 1.85 /mi feature has been examined 
in individual spectra of this source contained in the LRS data- 
base. None of the scans show anomalies at this wavelength or 
cosmic ray spikes. It is significant that the temperature of the 
blackbody that estimates the continuum of 15084 — 5701 is 410 
K, just the threshold temperature for hydrogen effusion in the 



Fig. 4 — xhe average normalized emission excess spectrum of five carbon 
stars 08189 - 3602, 15084 - 5702, 18320 - 0352, 18410 + 3654, 18551 +0323. 
These stars appear to have narrow emission features in the 11.3 to 12.0 ^m 
region that are associated with unsaturated (sp 2 ) bonds. 

laboratory. When this star changes into a planetary nebula, 
then the 11.3 /xm band should become apparent. We propose 
this circumstellar grain evolution is analogous to the labo- 
ratory experiments on a :C — H by Dischler et al. (1983a, b). 

Three other stars that are listed in Table 2 are interesting to 
note in regard to possible narrow PAH features: HK Lyr = 
18410 + 3654, 18551+0323, and 15094-6953. 18320-0352 
may be interesting but has low S/N. Their spectra are also 
shown in Figure 4. These spectra display peculiarly long 
central wavelengths of a narrow emission peak superimposed 
on the wide petestal. The narrow feature is not attributable to 
SiC or the typical PAH/UIR emission. Its width and center 
wavelength are variable form source to source. Both param- 
eters fall within an acceptable range for interpretation similar 
to that discussed for 15084-5701 in the preceeding paragraph. 
While there is an atomic line at 1 1.77 /xm due to Cl iv, it is not 
expected to be present in these objects because it requires a 
highly ionizing radiation. That radiation would produce other 
emission lines normally seen in planetary nebulae, that are not 
seen in these objects. Some of these objects belong to the 
coldest classes of carbon stars, and yet only a few objects from 
these classes display such features. Apparently low tem- 
peratures are required for the formation of PAH molecules, as 
well as some other unique circumstances that are as yet unde- 
termined. 


TABLE 2 

Peculiar Narrow 1 1 + ^m Bands in Carbon Stars 


LRS LRS 11+ /im 

IRAS Name Other Name ACT! AutoClass Class 2 ceBl<r 


15084-5701 ... 5 220:1 42 11.85 

15094-6953 X TrA 3 225 42 11.73 

18320-0352 RAFGL 7012S 5 220:1 42 11.85 

18410+3654 HK Lyr 3 225 42 11.43 

18551 +0323 RAFGL 2287 5 220:1 42 11.85 
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The 1 1.85 pm transition of PAHs is due to the C— H out-of- 
plane deformation mode of saturated edge sites. With an suffi- 
ciently cold carbon star, the 1 1.85 pm band ought to be appar- 
ent if these ring molecules are present and abundant. Its 
appearance at a temperatures less than 4000 K would indicate 
that oc:C— H edge sites become fully saturated below this tem- 
perature. Heating studies of a :C~ H samples in the laboratory 
that are made from radio frequency plasma deposited benzene 
on fused quartz and germanium substrates demonstrate the 
evolution of the unsaturated edge sites ( sp 3 ) to saturated sites 
(sp 2 either olefinic or aromatic; Dischler et al. 1983a, b) that is 
attributable to hydrogen effusion at temperatures above 400 K. 
The concentrations of these sites are modified by thermal 
annealing in two separate temperature regimes: 450-600 K 
and 600-900 K. Up to 600 K, the sp 3 /sp 2 ratio increases indi- 
cating the removal of sp 3 unsaturated bonds by hydrogen 
attachment and cross-linking. In the 600-900 K region, hydro- 
gen is effused until at 900 It, all the carbon atoms have sp 2 
aromatic hybridization. Such hydrogen effusion would mani- 
fest itself spectrally through an evolution of C— H emission 
features in an appropriately chosen sequence of objects. In the 
11 pm region, the sequence is ordered from ~~ \2 pm into 11.3 
pm features. The annealing process is not reversible in tem- 
perature, and implies that the starting material sp 3 
(unsaturated bonds) is created at low temperatures and subse- 
quently becomes sp 2 (saturated bonds). This sequence is in line 
with the supposed evolution of carbon stars into planetary 
nebular. It is not completely analogous to the induced nucle- 
ation theory proposed by Frenklach et al. (1989). 

We can conclude that when PAH or HACs are present in 
classical carbon stars, and they emit a narrow feature in the 1 1 
pm region, it occurs at about 1 1.9 pm, not 11.3 pm. In environ- 
ments where copious UV radiation is present, they emit at 1 1.3 
pm, not 11.9 pm. The evolution of the feature from 1 1.9 to 1 1.3 
pm reflects the effusion of hydrogen from the starting material, 
and the evolution is a result of thermal processing of the 
material. Material with unsaturated bonds evolves into that 
with saturated bonds, either from thermal or UV radiation 
induced heating. It is not necessary for PAHs to be generated 
in carbon star atmospheres and supplied directly to the inter- 
stellar medium (ISM). HACs in the ISM can be modified under 
the influences of UV radiation when their temperatures are 
elevated over the threshold of hydrogen effusion, 600 K, by 
absorption of UV photons. 

Pure carbon dust with C — C and C=*C bonds was not 
found. Graphite lacks identifiable bands in the LRS range. C 60 
displays two prominent bands at 17.55 and 19.3 pm (Onari, 
Tada, & Toshihiro 1991) that are not observed in any of the 
Class spectra or the individual object spectra. 

Below 500 K, an emission at k > 20 pm is noticeable in the 
spectra of Classes 1 and 5. The feature is the blue end of the 30 
pm feature that has been observed in carbon stars and planet- 
ary nebulae (Forrest et al. 1981 ; Goebel & Moseley 1985) and 
has been associated with MgS (Goebel 1980; Goebel & 
Moseley 1985; Nuth et al. 1985). This rise in emission is not to 
be confused with the 21 pm emission (Kwok, Volk, & Hrivnak 
1989) found in certain carbon stars with bright visual com- 
ponents, and associated with SiS 2 by Goebel (1993). This emis- 
sion at A > 20 pm is noticeable in the LRS spectra of individual 
stars, e.g., IRC + 10216, in the Classes in which the MgS band 
has already been observed and is not an artifact of detector 
saturation, data processing problems, or erroneous cali- 
brations. In Classes 1 and 5 we find an extensive observing list 


of stars that should be valuable to future airborne or satellite 
observing programs in the far infrared. 

The broadband colors of Classes 1 and 5 display the 
expected characteristics of the MgS band. In the [12]— [25], 
[25]— [60] color-color diagram, the Class mean color is shifted 
to the right, and in the [25]-[60], [60]— [100] diagram the 
Class mean color is negative . This negative color can be inter- 
preted as a color deficiency at 100 pm, i.e., as an emission 
feature that shows an excess in the 25 and 60 pm bands, and 
which returns to the continuum at 100 pm. The mean colors of 
Classes 1 and 5 are strikingly different than those of Class 0, 
which includes much hotter stars, and Class 6, which contains 
hot stars with significantly positive [60]-[100] colors 
(emission excesses) that can be interpreted as cold shells of 
carbon grains. Therefore, the colors of Classes 1 and 5 are not 
consistent with cool circumstellar shells of pure carbon grains. 

3.2. Evolution 

The observed correlation of Class temperature with the pres- 
ence of the SiC and a:C— H dust species suggests a grain 
condensation sequence in which SiC, the most refractory of 
carbon rich dust species, condenses out of the atmosphere first. 
As the atmosphere evolves, cools and expands, SiC forms the 
condensation nuclei for the subsequent formation of a:C— H 
(Martin & Rogers 1987). At still lower temperatures and later 
times, MgS catalysis occurs on the grains (Goebel & Moseley 
1985). Although the evolution of carbon stars is not clearly 
understood, this condensation sequence should occur in 
appropriate carbon stars whose circumstellar shells evolve 
from higher to lower temperatures. From the statistics of 
group membership in the Classes included in the sequence 
shown in Figure 1, roughly j (60%) of all the carbon stars that 
are included in this study can be identified as falling into this 
evolutionary sequence, There are a substantial number of 
spectra that do not fit into this sequence and display molecular 
absorption features, specifically, Class 3. 

It is significant that, as carbon stars are thought to be the 
red giant evolutionary phase that precedes the protoplanetary 
and planetary nebula phases, the spectral signatures of 
a:C“H are found at an earlier evolutionary phase than are the 
spectral signatures of PAHs. Although not included in the 
carbon stars studied here, the prototypical carbon rich proto- 
planetary nebula 04395 + 3601 = PK 166-6.1 = RAFGL 618 
(MetaClass yO) displays the broad absorption features without 
the narrow 1 1.3 pm emission, and yet has weak narrow PAH 
type features at 3.3 pm. The presence of the narrow emission 
features at 3.3 pm may be the result of nacient bipolar nebula 
formation in this object, or they may be a:C“H emission 
similar to that observed by Dischler et al. (1983a, b) in the 
laboratory. Since the dust temperature of RAFGL 618 is very 
low, ^150 K, it is unlikely that the PAHs exist as free mol- 
ecules, but rather have condensed completely into the solid 
phase. The vapor pressure of a heavy molecule like a PAH at 
150 K is absolutely negligible, and there is no evidence that 
PAHs sublimate. So before the planetary nebulae phase begins, 
the balance in the population of species should be dramatically 
in favor of simple C— O molecules, and not PAHs. As the 
planetary nebula develops, excitation of PAH emissions 
should progress and eventually dominate the spectrum. As an 
aside remark, Jourdain de Muizon et al. (1990) suggest that 
CRL 618 displays a silicate emission spectrum in the LRS 
range of wavelengths. We interpret the spectrum as displaying 
absorption bands of the molecules HCN and C 2 H 2 centered at 
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7 and 13 /an. CRL 618 is engaged in form of spectral mimicry, 
of which more will be said later. It is not the sole example of 
this phenomenon; many such objects are found in Class 3. 

Confirmation for this line of reasoning is found by consider- 
ing the spectra of the LRS 92 objects 08189-3602 = PK 
245 + 0.1, and 12063-6259 = He 2-77 = PK 298-0.1 = 
RAFGL 4144 (Fig. 5), members of the yO MetaClass, and of the 
same split Class y0:2 as RAFGL 618. All have a similar T c ~ 
180 K, but 08189-3602 and 12063-6259 have a weak 11.6 
pm narrow features on a broad 1 1 + pm pedestal along with 
strong 8.5 and 7.7 /mi emissions, as well as emission lines at 
12.9, 15.5, and 18.8 /on, due to, respectively, Ne n, Ne in, and 
S hi (Jourdain de Muizon et al. 1990). The latter two objects 
are probably not planetaries, rather they are most likely H ii 
regions (Simpson & Rubin 1990). They are too luminous by 1 
or 2 orders of magnitude based on their radial velocities and 
Galactic latitudes. They are relevant to the present discussion 
by virtue of their displaying unusually long central wave- 
lengths for their narrow 1 1 + /an emissions. Probably, PK 
245+0.1 and 298 — 0.1 have just started to dissipate their 
shells, and the 11.6 /an band should in time shift into an 11.3 
/an band. The dissipation process for protoplanetaries and 
young H n regions is probably very similar, with the details 
being controlled by the luminosity of the central source and its 
temperature. The lower limit to protoplanetary shell tem- 
peratures of carbon-rich objects prior to the onset of the plan- 
etary nebula eruption, e.g., emergence of the bipolar emitting 
regions that are presumably in the polar regions of the circum- 
stellar shell, is - 150 K. RAFGL 618 could exemplify the cir- 
cumstellar shell just prior to eruption. Kwok & Feldman 
(1978) observed RAFGL 618 in free-free emission at 6 cm and 
classified its central source as B2. The transition spectral type 
may be between B2 and 07 (A. G. G. M. Tielens 1994, private 
communication). See Martin-Pintado et al. (1992) for a sche- 
matic of RAFGL 618 deduced for NH 3 observations that is in 
line with this interpretation. 

Do our results imply an evolutionary scheme significantly 
different from previous proposals, or is an alternative scheme 


required? We are able to fit a large fraction of the carbon stars, 
those with circumstellar emission features, into a continuum 
temperature sequence. What we can say is this: // carbon stars 
evolve to lower temperatures and larger C/O, then the sequence 
of LRS spectral types in Figure 1 can be interpreted to represent 
an evolution in the circumstellar dust shell from a-SiC to a:C— H 
that makes physical sense. Due to the large number of free 
parameters, there are always exceptions to the rules. The 
crucial parameter L* is not well-determined for carbon stars of 
any class. Until their distances are determined reliably, L* and 
M* influences on the sequence in temperature will remain an 
open question. Because of this uncertainty, our simple 
sequence in temperature remains open to interpretation, yet 
the sequence remains established. 

We note that the LRS spectra alone are insufficient to disen- 
tangle the basic stellar parameters T eff , C/O, Af*, L*, etc. 
There have been recent studies examining the masses of cir- 
cumstellar envelopes of carbon stars and distances derived. 
The study of Olofsson et al. (1993) summarizes the available 
data for some of the stars (~ 120) in the present study. In order 
to arrive at distances, an absolute magnitude in the K band of 
— 8.1 was assumed. If all the stars have the same absolute 
luminosity, then the sequence in temperature is easily con- 
verted into evolution. However, that assumption of similarity 
is suspect when applied to our larger sample. 

Chan and Kwok’s evolutionary scheme can be stated simply 
as follows: Most carbon stars are formed before they begin to 
form dust shells. Thereafter , they evolve to lower temperatures 
and thicker SiC dust shells. The driving mechanism here is an 
increasing C/O that is the result of helium shell flashes in the 
interior and plumes that mix the triple-oc generated material to 
the surface, where it is driven outward into the circumstellar 
regions by a stellar wind. Chan and Kwok’s conclusions 
regarding the evolution of carbon stars is based largely on their 
comparison of radiative transfer models with the photometry 
in the literature that ranges from near to far infrared, a tech- 
nique similar to that employed previously by Rowan- 
Robinson & Harris (1982, 1983), van der Veen & Habing 




Fig. 5 a Fig - 56 

Fig. 5.— Spectra of the H n regions 08189-3602 ( a ) and 12063-6259 (6). These objects appear to have narrow emission features at 1 1.6 fim, rather than the 
nominal value of 1 1.3 fim associated with the UIR/PAH series of features. 
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(1988), and Willems & de Jong (1986a, b). The latter reached 
somewhat different conclusions regarding the evolution of 
carbon stars while employing similar data sets and models. 

Willems and Jong’s scenario is this: Most carbon stars are 
formed after dust shell formation. A further increase in C/O 
produces an evolutionary sequence in dust population in the shell 
starting at silicates and progressing to SiC and carbon. The 
scheme of Willems and de Jong recognizes the population of 
Class 1 1 and wrestles with the consequences. Chan and Kwok 
(as well as Zuckerman 1993) dismiss Class ll’s characteristics 
as extraneous or spurious. This difference largely accounts for 
the orthogonal viewpoints each group achieves with essentially 
the same starting data and analysis approach. 

Here, the evolutionary scenario for the dust outlined differs 
significantly from that proposed by Chan & Kwok (1990). 
They studied the circumstellar emission spectra of many of the 
same objects and concluded that the evolution of SiC from 
hexagonal-rhombohedric form, a-SiC, to the cubic form, 
ft- SiC, occurs at a transition temperature ~ 2000 K, and that 
the 8.7 pm emission feature is due to compound carbon molec- 
ular bands near the photosphere. They chose to describe the 
dust opacity throughout the infrared as “empirical SiC” in 
spite of the large body of research that had clearly shown that 
separate carbon and SiC populations are the appropriate 
description of the circumstellar dust. This distinction is more 
than mere semantics. The distinction between the two grain 
populations leads to different possibilities for describing not 
only the spectra, but also the evolutionary history of the shells 
and stars within them. The empirical SiC dust is distinctly 
different from the Si polluted carbon dust discussed by Martin 
& Rogers (1987), in that the latter considered the carbon dust 
to have occasional Si“C bonds incorporated in the bulk that 
gave rise to the SiC feature, while the C~ C bond complexes, 
whatever their exact form, produced a featureless continuum 
throughout the infrared. Carbon dust is well-founded theoreti- 
cally and observationally, especially for the infrared bright 
carbon star IRC +10216 (see Gail & Sedlmayr 1987). Chan 
and Kwok would have been better served by naming their 
opacity empirical carbon star dust rather than empirical SiC. 
The recent article by Lorenz- Martins & Lefevre (1993) thor- 
oughly discusses the multicomponent dust model through 
radiative transfer calculations and compares the results with 
observations. The weight of evidence is clearly on the side of 
multiple dust components. 

Our Class 11, those stars that display silicate like emission, 
or its equivalent in other studies, is at the center of the debate 
in the literature relating to the evolutionary paths of carbon- 
rich stars and their circumstellar shells (Zuckerman & Madda- 
lena 1989; Zuckerman 1993; Willems & de Jong 1988; de Jong 
1989; van der Veen & Habing 1988). The resolution of the 
debate must explain Class 11, not ignore it. It is outside the 
scope of this work to address this issue beyond pointing out 
that Class 11 will not go away and hide; it demands explana- 
tion. We do point out the possible role of spectral mimicry 
apparent in some carbon star spectra, such as T Lyr and 
AFGL 618. The spectra of some, but not all, of the objects in 
Class 1 1 are attributable to this mechanism. 

We suggest that our spectral sequence based on T c mimics 
an evolutionary sequence in time and/or mass paralleling an 
increasing C/O. If the evolution of carbon stars from visually 
bright to infrared bright is indeed correct, then the evolution 
must be toward increasing C/O. We consider the set of LRS 
spectra to reflect evolution in time of carbon stars and to rep- 


resent the various states they would possess if sampled ran- 
domly. Then we can consider the progression of spectral types 
in Figure 1 to recapitulate the evolution of such stars. This 
interpretation is analogous to the tranditional spectral classi- 
fication sequence in visible astronomy that uses a two- 
dimensional parameter set consisting of temperature and 
spectral feature strength. We propose that such a two-dimen- 
sional parameterization also applies to dust features in the 
infrared spectra of carbon stars and is a useful classification 
methodology. The two parameters are T c and feature strength 
of the SiC and a:C — H components. Because broadband 
color-color diagrams are not very sensitive to the feature differ- 
ences, the use of these diagrams does not produce a unique 
mapping into these two parameters, so it is necessary to utilize 
spectrophotometry as is present in the LRS. A narrowband 
filter set specifically tuned to the dust features would give 
useful information. The selection of such a set would require 
careful consideration due to the presence of Class 3 stars in any 
random sample of carbon stars. 

3.3. Dust Condensation and Mass Loss 

Homogeneous dust condensation theories require carbon 
dust to nucleate on preexisting molecules and cluster complex- 
es, not on condensation nuclei like SiC or TiC. From homoge- 
neous condensation calculations, it has become well known 
that SiC condenses at a higher temperature than carbon dust 
(Gilman 1969; McCabe 1982). Class 4 is consistent with homo- 
geneous condensation of SiC. Various types of carbon dust 
(graphite, amorphous carbon, tetrahedrally coordinated 
carbon, hydrogenated amorphous carbon, diamonds, etc.) are 
considered to form depending upon a variety of initial condi- 
tions. There are identifiable instances where homogeneous 
nucleation of carbon dust is consistent with the observations, 
such as the hydrogen deficient R CrB stars. Because the chem- 
istry of hydrocarbons is very complex, there have been a wide 
range of proposals as to the precise molecular pathway (C 3 , 
C 2 H 2 , PAH, etc.) from which the dust (a:C— H, graphite, 
HAC, amorphous or glassy carbon) forms. The red giant phase 
should be warm enough (1500 < T e < 3000 K) to allow diverse 
populations of molecules and dust as has been modeled by 
Gail & Sedlmayr (1987). Their models include PAHs in a ther- 
modynamically natural manner. For most carbon stars in our 
study, we have found that homogeneous dust populations are 
not the rule, and heterogeneous populations predominate. 

From other studies there is mounting evidence that dust 
condensation in carbon rich circumstellar shells proceeds het- 
erogeneously. Studies by McCabe (1982), using equilibrium 
nucleation theory, Gail & Sedlmayr (1984) and Frenklach et al. 
(1989) conclude that SiC condenses out closer to the photo- 
sphere and at higher temperatures than does amorphous 
carbon dust. Hence, as mass is driven outward from the photo- 
sphere deeper into the circumstellar shell by the stellar wind, 
where lower temperatures prevail, SiC is injected as readily 
available nuclei or seedlings for the heterogenous condensation 
of carbon. This scenario is consistent with our interpretation of 
the class spectra and needs theoretical confirmation in regard 
to the essential details. While it is clear that there should be 
spectral changes due to the incorporation of various species, 
one might suspect structural differences in the circumstellar 
shells. In the models, the resulting global thermodynamic and 
hydrodynamic shell structure, appears not to be so much de- 
pendent on the particular dust species, but rather the fact that 
dust is present or absent (Gail & Sedlmayr 1987; Fleischer, 
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Gauger, & Sedlmayr 1992; Winters, Dominik, & Sedlmayr 
1994). 

Heterogenous condensation may have implications for the 
explanation of certain anomalous gas phase molecular abun- 
dances. Significant amounts of Si can be consumed by SiC. 
Observations of Goldhaber & Betz (1984) and Keady & 
Ridgway (1993) point to an anomalous underabundance of 
silane inside 40 R * that is unexplainable in terms of equi- 
librium or nonequilibrium gas phase chemistry in or near the 
photosphere. The Berlin group models show that condensation 
of carbon dust in carbon rich atmospheres usually proceeds 
through C 2 H 2 or C 2 H implying hydrogenated carbon of some 
form is a likely condensate (see also Keady & Hinkle 1988). 
HAC mantles on SiC grains could make the Si unavailable for 
gas phase or grain surface chemistry. This is an interesting area 
for further theoretical work. 

The recent study of Lorenz-Martins & Lefevre (1993) exam- 
ines the composition of a two component dust with radiative 
transfer models and compares the resultant flux curves with 
nine carbon stars. They find a range for the number density of 
grain populations, N SiC JN AC of 0.03 to 0.20, and support the 
conclusion of Skinner & Whitmore (1988a, b) that the strength 
of the SiC feature is related to the total mass-loss rate. Skinner 
& Whitmore (1988a, b) related the strength of the SiC emis- 
sion, P(SiC), expressed as the product of the band emission to 
continuum ratio, the IRAS band I emission intensity, and dis- 
tance squared of the star, directly to the CO radio emission line 
intensity and thereby the mass-loss rate, dM/dt . Lorenz- 
Martins & Lefevre (1993) related optical depth, t cx1 , and 
N SiC /N AC . All four parameters appear to be related to the mass 
loss rate in the circumstellar shell. The correlation between i ext , 
N SiC /N KC , P(SiC), and CO integrated line intensity appears to 
be positive, allowing for a few exceptional stars like RY Dra. 
Neither Skinner & Whitmore (1988a, b) or Lorenz-Martins & 
Lefevre (1993) recognized the difference between the SiC band 
and the ot:C— H band, nor that of the “windowing effect ” of 
molecular absorption bands seen against hot continua, a 
coincidence of absorption features that we refer to a spectral 
mimicry . The distinction does not invalidate their result; 
rather, it may account for some of the dispersion in their mea- 
sures. For example RY Dra is a member of carbon star Class 3, 
a Class that is dominated by hot continua with molecular 
bands in absorption. By shear coincidence, spectral mimicry , it 
happens that the RY Dra emission feature is the window 
between two HCN absorption bands (see Fig. 6). It is not really 
valid to place it in Figure 12 of Lorenz-Martins & Lefevre 
(1993) that correlates r ex , and N SjC /jV AC . Extraction of RY Dra 
improves the correlation between those two parameters. In 
Figure 6, the 1 to 25 pm infrared spectrum of RY Dra is shown 
to exemplify the spectral mimicry concept. For similar reasons, 
the exclusion of Y CVn from the study of Skinner & Whitmore 
(1988a, b) could improve the scatter in their relationships. In 
fact, all Class 3 objects should be carefully considered in such 
studies to avoid confusion by spectral mimicry. 

Although Lorenz-Martins & Lefevre (1993) did not use core 
mantle grains in their study, they favor an interpretation of the 
data in such terms, and promised a future contribution relating 
to it. Such heterogenous nucleation has been favored by 
several investigators of carbon star dust condensation (e.g., 
McCabe 1982; Frenklach et al. 1989). However, the theory is 
most commonly applied to the condensation of ice mantles in 
dark clouds (e.g., Tielens & Allamandola 1987), and the results 
are not directly applicable to the conditions in carbon star 
circumstellar shells. 



Fig. 6.— The normalized emission excess spectrum of the carbon star RY 
Dra between 1.2 and 24 /im. This figure combines spectra taken from NASA’s 
Kuiper Airborne Observatory and NASA’s IRAS mission. The spectra were 
obtained on the dates: 1. 2-5.6 /an, 1976 December 14; 4-8 //m, 1979 April 10; 
8-24 /an; 1983 January-November. Continuum removal included the contri- 
butions of the photosphere and a 300 K. amorphous carbon shell. 


Whatever their origins, the strengths 1 1 4- pm emission fea- 
tures appear to be related to the total mass-loss rate, as long as 
the features are optically thin. As the features express the abun- 
dances of different dust populations, the exact nature of the 
correlation is in need of revision. There are probably two 
tracks in Figure 1 of Skinner & Whitmore (1988a, b), one for 
pure SiC, corresponding to Class 4, and the other for pure 
a:C— H, corresponding to Class 5, as well as a continuum of 
mixed states, corresponding to Classes 0 and 1 . 

We anticipate that the SiC emitting stars, Class 4, should 
show a lower mass-loss rate than the a:C— H stars. Class 5, 
because the latter tend to be the infrared bright carbon stars, 
and the former the visually bright carbon stars. Indeed, there 
are no members of Class 4 in the sample of Skinner & Whit- 
more (1988), indicating a lack of observable CO radio emission 
lines from which distances were estimated. There are 10 
members of Class 0, i.e., those with mixed SiC and a:C— H and 
continua temperatures averaging 800 K. These 10 members of 
Class 0 that are common to both studies have an average of 
dM/dt ~ 4.2 x 10 -6 M q yr _1 . The values are dispersed over 2 
orders of magnitude. The two members of Class 5 are IRC 
+ 10216 and CRL 2155, both of which exceed dM/dt ~ 
2 x 10 -5 Mq yr' 1 . 

The Class 3 stars, also visually bright, should show measur- 
able mass-loss rates. The dust-to-gas mass must be measured 
by an alternate procedure that takes into account the strength 
of the continuum of an emitting dust shell. This must be esti- 
mated by subtraction of a photosphere, after careful consider- 
ation of the effects of molecular features. Indeed there are very 
few stars for which there is available sufficient spectro- 
photometric data spanning the 1 to 30 pm spectral region for 
the conduct of such a study. The current authors are preparing 
such a compilation using aircraft and satellite data. The spec- 
trum of RY Dra shown in Figure 6 is an example. The 
members of Class 3 in the sample of objects used by Skinner & 
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Whitmore 1988a, b) have consistently low dM/dt ~~ 1.9 x 10“ 7 
M © yr ” l . CIT 6, not included in this estimate, is atypical of the 
sample of seven sources and has a dM/dt ~ 3.2 x 10“ 6 Af© 
yr“ l . It may have classification problems induced into its spec- 
trum by the nonlinear detector effects of extremely bright 
sources. IRC -1-10216 clearly has this problem. 

Based upon this small sample of common objects in the two 
studies, it appears that the division of spectra and association 
of spectral features proposed here is reflected in the mass-loss 
rates of the objects as measured by Skinner & Whitmore 
(1988a, b). Generally, circumstellar shells with colder tem- 
peratures display a:C— H and have higher mass-loss rates; 
warmer circumstellar shells display SiC and have lower mass- 
loss rates. There are two notable exceptions, CIT 6 and IRC 
-I- 10216. Exceptions are to be expected in any study of carbon 
stars. Hopefully improvements in radio receivers will result in 
sufficient sensitivity for the measurement of distances necessary 
to this study. A further study of the entire sample of carbon 
stars in the LRS with these more sensitive radio receivers 
ought to clarify this division. 

The combination of homogeneous condensation of seedlings 
and heterogeneous condensation of mantles has been termed 
induced nucleation by Cadwell et al. (1994). Their model elabo- 
rates the proposition by Frenklach et al. (1989) that a circum- 
stellar shell surrounding a carbon-rich red giant star has 
multiple zones in which condensation takes place homoge- 
neously nearest the star, and heterogeneously for amorphous 
carbon and PAH condensates in more extended zones. The 
discussion of our spectral sequence lends credence to such a 
scheme by identifying Classes of carbon star spectra consistent 
with such zones. Class 4 has the hottest continuum tem- 
perature, 3000 K, and shows exclusively SiC emission. Class 0, 
the next hottest at 800 K, shows the beginnings of the 8 -f /*m 
a:C— H band and the SiC band combined with the 1 1 + /im 
a:C— H band. Classes 1 and 5 are progressively colder, 600 
and 470 K, and exhibit stronger a:C“H bands, and little to 
none of the SiC band. Finally, two objects, 15804 — 5701 and 
08189 — 3602, exhibit unsaturated (sp 3 ) hydrocarbon emission 
bands at 11.85 and 11.6 pm and very low continuum tem- 
peratures of 410 and 1 50 K. 

We could consider the spectra to sample progressively 
colder regions of the circumstellar shell as those regions 
approach optical depth one progressively farther from the 
photosphere. The temperatures noted in Figure 1 do not agree 
precisely with those quoted by Frenklach et al. (1989) in their 
schematic diagram of the circumstellar shell's dust forming 
regions. Our temperatures are continuum emission tem- 
peratures, and may not accurately reflect the conditions in the 
formation region for a particular species, but rather that of the 
dominant emitting region underneath. A detailed radiative 
transfer computation is required to associate radiation tem- 
peratures with features in the emitting spectra. In spite of this, 
the sequence of species deduced is in agreement with their 
schematic. 

Heterogenous nucleation, or induced nucleation, or catalysis 
is almost certainly required to explain the presence of MgS and 
SiS 2 emission features in many carbon-rich circumstellar shells 
(Goebel & Moseley 1985; Nuth et al. 1985; & Goebel 1993). 
Because of the relatively low abundances of the sulfides com- 
pared with carbon, nearly all the sulfur in the coldest regions of 
the shells must condense out onto relatively large grains. This 
is only possible on preexisting grains of carbon through, for 
example, an efficient catalytic reaction. It is worthwhile to note 


that although there are many instances of MgS band (30 /im) 
indicated in the long-wavelength spectra of the carbon stars in 
our sample, there are no occurrences of the more easily per- 
ceived SiS 2 band (21 /im). If the 21 /an band were due to the 
presence of some hypothetical carbon compound, we find its 
lack of occurrence difficult to explain. 

Heterogenous nucleation, or induced nucleation with SiC 
seedlings, is limited in its application to carbon stars. It cannot 
be invoked to explain the condensation process in hydrogen 
deficient R CrB type stars. They show no SiC bands, so the 
nucleation must be homogenous, at least in the sense that only 
carbon compounds are involved. These stars are thought to 
have much higher surface gravity and therefore, gas densities, a 
situation that enhances comparatively three-body collisions, a 
theoretical barrier to homogeneous nucleation theory. 

4. CONCLUSIONS 

A new classification scheme of the IRAS LRS carbon stars 
has been presented. It comprises the separation of 718 prob- 
able carbon stars into 12 distinct self-similar spectral group- 
ings. Continuum temperatures are assigned and range from 
470 K to 5000 K. Three distinct dust species are identifiable: 
SiC, a:C— H, and MgS. SiC and a:C— H band strengths are 
found to correlate with the temperature progression among 
Classes. If temperature corresponds to evolution, then SiC 
grains precede oc:C — H and form the nuclei for the conden- 
sation of the latter material. At still lower temperatures, MgS 
appears to be quite prevalent. We conclude that SiC seeding of 
the circumstellar gas is a prerequisite for the efficient growth of 
hydrocarbon grains. Heterogenous (or induced) nucleation 
appears to be a common, although not an universal, mecha- 
nism in carbon rich circumstellar shells. It is thermodynami- 
cally efficient and able to account for the multitude of species 
that give rise to the several 1 1 + fim spectral features apparent 
in these objects. 

No 11.3 /im PAH features are identified in any of the 718 
carbon stars. However, a handful of objects, IRAS 
15048 — 5702, display an 11.85 /im emission feature character- 
istic of laboratory samples of coronene, a PAH. That feature 
corresponds to the C— H vibrational waging of saturated edge 
sites. These few examples of narrow emissions between 11.3 
and 12 /an are noted and attributed to hydrocarbon bonds for 
the first time. The following puzzle presents itself. If a:C— H or 
HACs are developed from PAH molecular constituents, and 
these must be very numerous, then why are not the PAH 
narrow emission features seen in the spectra of any of the 
hundreds of carbon stars from the LRS Atlas? One UV photon 
corresponds energetically to about 100 each 1 1 /im photons. 
Certainly, the density of PAHs in a circumstellar environment, 
must be more than 100 times greater than that in the inter- 
stellar medium. oc:C— H condenses at higher temperatures 
than do the aromatic hydrocarbons (PAHs), and this fact may 
explain this situation. 

The new classification scheme yields a progression of spectra 
which when combined with estimated temperature of the 
underlying continua, can be interpreted as an evolutionary 
progression in time or mass. A two-dimensional classification 
of temperature and feature strength can be extracted that cor- 
responds to this evolution. It is vitally important to determine 
the distances to many carbon stars individually and reliably so 
that their luminosities can be determined. Only then can the 
influence of time and mass be determined on this evolutionary 
sequence. 
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Note added in proof . — A similar discussion regarding dehydrogenation of carbon analog grains has appeared in print 
(V. Mennella et al., ApJ, 444, 288 [1995]) in its application to the 2200 A band in the interstellar medium. 




